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In a celebrated paper, Dyson shows that the spectrum of an 
n X n random Hermitian matrix, diffusing according to an Ornstein- 
Uhlenbeck process, evolves as n noncoUiding Brownian motions held 
together by a drift term. The universal edge and bulk scalings for 
Hermitian random matrices, applied to the Dyson process, lead to 
the Airy and Sine processes. In particular, the Airy process is a con- 
tinuous stationary process, describing the motion of the outermost 
particle of the Dyson Brownian motion, when the number of parti- 
cles gets large, with space and time appropriately rescaled. 

In this paper, we answer a question posed by Kurt Johansson, 
to find a PDE for the joint distribution of the Airy process at two 
different times. Similarly we find a PDE satisfied by the joint dis- 
tribution of the Sine process. This hinges on finding a PDE for the 
joint distribution of the Dyson process, which itself is based on the 
joint probability of the eigenvalues for coupled Gaussian Hermitian 
matrices. The PDE for the Dyson process is then subjected to an 
asymptotic analysis, consistent with the edge and bulk rescalings. 
The PDEs enable one to compute the asymptotic behavior of the 
joint distribution and the correlation for these processes at differ- 
ent times ti and t^, when t2 — ti — > oo, as illustrated in this paper 
for the Airy process. This paper also contains a rigorous proof that 
the extended Hermite kernel, governing the joint probabilities for the 
Dyson process, converges to the extended Airy and Sine kernels after 
the appropriate rescalings. 
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2 M. ABLER AND P. VAN MOERBEKE 

1. Stating the results. The Dyson Brownian motion [4] 

(Ai(t),...,A„(t))GM" 
with transition density p{t, fi, A) satisfies the diffusion equation 

dt 2^dXi ^ ^dXi^xy 



^fl_^ d dlog^/W{X} \ 

VV2 9A? dXi dXi r 



with 



^{X) = ^\X)X{e 
1 



Roughly speaking, it represents n Brownian motions repelling one another, 
with the exponential in ^{X) having the effect of preventing the system 
from flying out to infinity. In his beautiful paper, Dyson generalizes the ran- 
dom matrix ensembles in such a way that the Coulomb gas model acquires 
a meaning as a dynamical system, rather than a static model. He shows the 
repelling Brownian motion above corresponds to the motion of the eigen- 
values (Ai(t), . . . ,Xn{t)) of a Hermitian matrix B, evolving according to the 
Ornstein-Uhlenbeck process 

I, J — i J 



with transition density (c : = e V^-^) 

P{t,B,B) = Z-'- L_e-( (1-^))) Tr(B-cBf 

(1 — c"^)" 

The -Bjj's in (1.1) denote the free (real) parameters in the Hermitian 
matrix B, with the Ba^s being its diagonal elements. In the limit t — > oo, 
this distribution tends to the stationary distribution 

n 

^-ig-(i/a^)TrB^ = Z-'A\X) n e-"?/'^' dX,. 

1 

With this invariant measure as initial condition, the joint distribution reads 

P{B{0)£dBi,B{t)edB2) 

(1.2) ^ r^^l dBidB2 (l/(a2(i-c2)))Tr(B?-2cBiB2+Bg) 

(l_c2)nV2^ 

for which a nonlinear PDE will be found in Theorem 1.1. According to 
[9], the joint probabilities for the Aj's can also be expressed in terms of a 
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Fredholm determinant of the so-called extended Hermite kernel Kl^l^{x,y), 
a matrix kernel, defined in Section 7. As elaborated in that section, for Ei 
and E2 C M, we have, for 1 <i <n, 

(1.3) P(an Xiih) G El, all A^ta) G E2) = det(I - (x^ci^£>^c)i<,,,<2). 

k i 

Since expression (1.2) is symmetric in Bi and B2, the probability (1.3) for 
the Dyson process is symmetric in Ei and £'2- Throughout the paper, we 
normalize the problem, by setting a= 1. 

The Airy process is defined by an appropriate rescaling of the largest 
eigenvalue A„ in the Dyson diffusion, 

(1.4) A(t) = lim V2n^/^(Xn(n-^^h) - V2^), 

n — >oo 

in the sense of convergence of distributions for a finite number of t's. This 
process was introduced by Prahofer and Spohn [12] in the context of poly- 
nuclear growth models and further investigated by Johansson [8]. Prahofer 
and Spohn showed the Airy process is a stationary process with continuous 
sample paths; thus the probability P{A{t) < u) is independent of t, and is 
given by the Tracy-Widom distribution [14], 

(1.5) P{A{t)<u) = F2{u):=exp(^- J^^ {a - u)q'^{a) da^ , 



with q{a) the solution of the Painleve II equation 

r g-{2/3)a3/2 

(1.6) q" = aq + 2q' with q{a) - "^7^^' « ^ 

[ y^—a/2, for a \ —00. 

Here, the joint probabilities for the process A{t) can also be expressed in 
terms of the Fredholm determinant of the extended Airy kernel K^f- {x,y), 
another matrix kernel, which is an appropriate limit of the extended Hermite 
kernel above; see [6, 8, 11, 12]. It leads to 

(1.7) P{A{h) G EuA{t2) e E2) = det(I - (XscK4xec)i<,,,<2), 

which is also symmetric in Ei and E2, as a consequence of the symmetry 
for the Dyson process. 

At MSRI (Sept. 2002), Kurt Johansson, whom we thank for introducing 
us to the Airy process, posed the question whether a PDF can be found 
for the joint probability of this process; see [8]. The present paper answers 
this question (Theorem 1.2), which enables us to derive the asymptotics of 
the large time correlations for the Airy process (Theorem 1.6); this question 
was posed in [12]. Our results on the Airy process for the special case of 
semi- infinite intervals appeared in [2], as well as the asymptotics. 
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The Sine process is an infinite collection of noncolliding processes Si{t), 
obtained by rescaling the bulk of the Dyson process in the same way as the 
bulk of the spectrum of a large Gaussian random matrix; namely, 

(1.8) Si{t):= lim K/2+i[ t;— ] for -oo<i<cx) 

in the sense of convergence of distributions for a finite number of t's. This 
process was defined by Tracy and Widom in [17]. Similarly, by taking the 
bulk scaling limit of the extended Hermite kernel (1.3), the joint probabilities 
for the SiS can also be expressed in terms of the Fredholm determinant of 
the extended sine kernel Kf,^,(x,y), yet another matrix kernel, 

(1.9) P(all S,(ti) G E'i, all S^it^) G i?|) = det(/ - ix,^K,%,X^J,^,^,^,), 

where here Ei and E2 must be compact for it to make sense. Note this 
probability is, as usual, symmetric in Ei and £'2- These kernels will be 
discussed in Section 7. For this process. Theorem 1.4 gives a PDF for the 
joint probabilities. 

The disjoint union of intervals 

r s 

Ei:=[J[a2i-i,a2i] and -E2 := IJ [^2i-i, M ^ 

1=1 i=l 

and t = t2 — ti specify linear operators, setting c = e~*, 

2r n 2s n 

2r ^ 2s Q 

2r Q 2s Q Q 

The duality ^ bj reflects itself in the duality Bi. We now state: 

Theorem 1.1 (Dyson process). Given ti < t2 and t = t2 — ti, the loga- 
rithm of the joint distribution for the Dyson Brownian motion (Ai(t), . . . , A„,(t)), 

(1.11) G„(t;ai,...,a2r;6i,...,62s) := log P(a// Xi{ti) G Ei, all Xi{t2) G E2), 

satisfies a third- order nonlinear PDE in the boundary points of Ei and E2 
and t, which takes on the simple form, setting c = e^*, 

(1.12) f^;^^^" =B^ ■^'^'^'^ 



(1.10) 



'BiAiGn + 2nc AiBiGn + 2nc' 
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The proof of this theorem will be given in Section 3. 
Similarly, the disjoint union of intervals 

r s 

(1.13) El := [J[u2i-i,U2i] and E2 := [J[v2i-i,V2i] ^R, 

i=l i=l 

and t = t2 — ti define another set of linear operators 

a ^ ^ d 



" ^ dui ' ^ ^ dvi ' 

We now give the equations for the joint probabilities of the Airy and Sine 
processes, which will be proven in Section 4: 

Theorem 1.2 (Airy process). Given ti < t2 and t = t2 — ti, the joint 
distribution for the Airy process A{t), 

G{t-Ui,...,U2r;Vi,...,V2s) := log P{A{ti) £ Ei, A{t2) £ E2) , 

satisfies a third- order nonlinear PDE in the Ui, Vi and t, in terms of the 
Wronskian {f{y),g{y)]y ■=f'{y)g{y) - f{y)g'{y), 

^-^ -j^g-j {{Lu + Ly)(LuE^ — LjjEu) + t^{Lu — Ly)LuL^)G 

= H(^l~ ^l)^^(^^ + ^^)'^^}L^+Lv■ 
COROLLARY 1.3. In the case of semi- in finite intervals Ei and E2, the 
PDE for the Airy joint probability, 



H{t;x,y) ■.= log p(A{ti) < ^,A(t2) < 



y-x 



V ' ' - 2 ' ' - 2 

takes on the following simple form in x,y and t^ , with t = t2 — ti, also in 
terms of the Wronskian, 

fl 16) 2i -(t^^ x^\(^'^^ d'^H\ Jd'^H d'^H') 

dtdxdy \ dx dyjydx^ dy"^ J [dxdy^ dy^ )y 



Remark. Note for the solution H{t;x,y), 



lim H{t; x, y) = log F2 ^min 



y + x y — X 



Theorem 1.4 (Sine process). For ti <t2, and compact Ei and E2 C M, 
the log of the joint probability for the sine processes Si{t), 

G{t;ui,...,U2r;vi,...,V2s) :=logP(a// Si{ti) £ E'i_, all Si{t2) £ E^) , 
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satisfies the third-order nonlinear PDE, 



(1.17) 



'u 



{2EyLu + {El, — Eu — l)L^)G 
(L, + L,)2G + 7r2 
_ J (2-EttL^ + {Eu — Ey — l)Lu)G 

(L, + L„)2G + 7r2 



Corollary 1.5. In the case of a single interval, the logarithm of the 
joint probability for the Sine process, 



H{t;x,y) = logP{S{ti) ^ [xi +X2,xi -X2],S{t2) i [yi +2/2,2/1 -2/2]) 



In a recent paper, Tracy and Widom [16] express the joint distribution of 
the Airy process for several times ti, . . . , t^j in terms of an augmented system 
of auxihary variables, which satisfy an implicit closed system of nonlinear 
PDEs. In [17], Tracy and Widom define the Sine process and find an implicit 
PDE for this process, with methods similar to the one used by them in 
the Airy process. The quantities involved are entirely different from ours 
and their methods are functional-theoretical; it remains unclear what the 
connection is between their results and ours. 

The PDEs obtained above provide a very handy tool to compute large 
time asymptotics for these different processes, with the disadvantage that 
one usually needs an assumption concerning the interchange of sums and 
limits; see Section 6 and the Appendix. This is now illustrated for the Airy 
process, for which we prove the following theorem, assuming some conjec- 
ture, mentioned below. This will be discussed in Section 6. A rigorous proof 
of expansion (1.19) was given later by Widom [18]; his proof was based on 
the Predholm determinant expression for the joint distribution. 

Theorem 1.6 (Large time asymptotics for the Airy process). For large 
t = t2 — t\, the joint probability admits the asymptotic series 



satisfies 



(1.18) 



d {2Ey d/dxi + {Ey -E^-l) d/dyi)H 

dxi {d/dxi+d/dyi)^H + Tr^ 

d {2E^d/dyi + {E,-Ey-l)d/dxi)H 

dyi {d/dxi + d/dyiYH + 7r2 



P{A{ti)<u,A{t2)<v) 



(1.19) 
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with the function q = q{a) as is the function (1.6) and 

(i 



^{u,v) := F2{u)F2{v) 



/ da\ i q^ da 
+ q\u)[\q\v)-\^j\^da 

poo rco 

+ / da{2{v - a)q^ + g'^ - q^) / q 



da 



Moreover, the covariance for large t = t2 — ti behaves as 
(1.20) EiAit2)Ait,)) - E{A{t2))E{A{t,)) = 1 + ^ + . ■ 
where 

c := 2 / / <t{u,v)dudv. 



Conjecture. The Airy process satisfies the nonexplosion condition for 
fixed x: 

(1.21) Ymi P{A{t)>x + z\AiQ)<-z) = {). 

z — ^oo 

This conjecture will be discussed in Section 6, just before the proof of 
Theorem 1.6 and in the Appendix. 

Finally, in Section 7, we give a rigorous proof of the convergence of the 
extended Hermite kernel to the Airy and Sine kernels, under the substitu- 
tions 

n 

, X I— > v2?T- + 1 H ^ — — 

S,:={t^'s^^ ^^^'^ 



n 



^/2?ii/6 



S2 ith^ ,s I— > ,x I— > ,y 



2n 2n ' ^2^'^ ^2^ j ' 

The precise formula for these kernels will be given later in the beginning of 
Section 7. 

Proposition 1.7. Under the substitutions Si andS2, the extended Her- 
mite kernel tends, respectively, to the extended Airy and Sine kernel, when 
n — > 00, uniformly for u,v £ compact subsets C M; 



hm K^f\x,y)dy, = e'^^" /'^^'"'^^ KfJu,v) dv. 
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2. The spectrum of coupled random matrices. Consider a product en- 
semble (Ml, M2) e nl := Tin X Tin of n X n Hermitian matrices, equipped 
with a Gaussian probabihty measure, 

(2.1) C„ dMi dM^e" Tr(M2+M|-2cMiM2)/2^ 

where dMi dM2 is the Haar measure on the product TC'^ , with each dMi , 

n n 

(2.2) dMi = Al{x)Y[dxidUi and dMa = A^(y) ]^ % dC/2 

1 1 

decomposed into radial and angular parts. In [1], we defined diff'erential 
operators Ak, of "weight" k, which form a closed Lie algebra, in terms 
of the coupling constant c, appearing in (2.1), and the boundary of the set 

r s 

(2.3) E = ElxE2■.= \J[a2^-l,a2^] x |J [62^-1, M C M^. 

i=l i=l 

Here we only need the first few ones: 



(2.4) 



2s Q Q 

In [1], we proved the following theorem in terms of the Wronskian {/, g}x = 
g{Xf) — f{Xg), with regard to the first-order differential operator X and 
based on integrable and Virasoro theory: 

Theorem 2.1. Given the joint distribution 

(2.5) Pn{E) := P{all (Mi-eigenvalues) G Ei, all {M2- eigenvalues) G E2), 

the function Fn{c;ai, ... ,a2r,bi, ... ,b2s) ■=logPn{E) satisfies the nonlinear 
third- order partial differential equation 

nc 



(2.6) 



B2AlFn,BiAlFn + 



C2-1 



nc 



A2BiFn,AlBiFn + ^^} =0 



C2-1 



Remark. Note that both Pn(£'i x E2) and PniE^ x E^) satisfy the same 
equation. 



PDES FOR DYSON AND AIRY PROCESSES 



9 



3. The joint distribution for the Dyson Brownian motion. 

Proof of Theorem 1.1. Using the notation of Section 2 and the 
change of variables 

(3.1) Mi= , ' ^=, 

one computes for t = t2 — > 0, with e~* = c, 

Gn{t] 0,1,..., a2r', bl,..., b2s) 

:= logP(an Ai(ti) G Si,an Xi{t2) G E2) 

= logP(aIl (i?(ti)-eigenvalues) G Ei,al\ (i?(t2)-eigenvalues) G E2) 
= \og jj ^-l_^:^l^^g-(l/(l-c2))Tr(B2+B2-2cBiB2) 



-eigenvalues S-Ei (1 — C^)" 

all _B2 -eigenvalues G_E2 



log / / 

J J a. 



all A/i -eigenvalues G-E'i/(^(l-c2)/2 ) 
all M2-eigenvalues ^E2/ {yj {l~c^)/2) 



Z'-^ dMi dM2 



Fn[c; 



^ g-Tr(M2+Af|-2cMiM2)/2 

ai a2r 



V(l-c2)/2''" V(l-c2)/2 
bl b2s 



\/(l-c2)/2''" V(l-c2)/2. 
in terms of the function Fn, defined in Theorem 2.1. Setting 

Fn{c;ai, . . . ,a2r]bi, . . . ,h2s) 

= Gn{t; aW{l-^)/2, . . . ,a2.V(l - c2)/2; 



bW{l-c^)/2, 62sV(l-c2)/2) 

= G„(t;ai, . . . ,a2r;bi, . . . ,62s) 

in (2.6) leads to the following equation for G„ := Gn{t; ai,. . . , a2r', bi,. . . , b2s), 
namely, 

BAGnAAGn + ^?^}, -UAGn,MlGn + ^?^}. =0, 



where 

'4.iF„(c; Qi, . . . , a2r', bi,..., b2s) =AiGn{t; ai, . . . , d2r', bi,. .. ,62s), 

BiFn{c; a2r; fel, • • • , &2s) =SiGn{t; Ol, • • • , a2r; ^1, • • • ,&2s), 
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with the Ai,Bi as in (2.4), and 

1 1 d 

Then clearing the denominators in (3.2) leads to 

(3.3) {B2AiGn,J3iAiGn + 2nc}^, = {^2^1 A^iCn + 2nc}B, , 

for the operators (1.10), with G.„ = G„(t; ai, . . . , 02r; &i, • • • , &2s)- This estab- 
lishes Theorem 1.1. □ 

Remark. In view of the remark in Section 2, here the expressions 

r logP(all Ai(ti) G ^i,an \(t2) G ^2) 

(3.4) G„ = <^ or 

[ logP(all Ai(ti) G ^f,all Ai(t2) G ^f) 

also satisfy the same equation. 

4. The joint distribution for the Airy process. 

Proof of Theorem 1.2. Consider as in (1.13), the disjoint union of 
intervals 

r s 

(4.1) El ■.= [j[u2i-i,U2i] and E2 := [j[v2i-i,V2i] <^^. 

i=l 1=1 

Then, setting t = T2 — ti, 

Hn{T;Ui, . . .,U2r;Vi, . . . ,V2s) 

_ / V2nV6(A„(n-V3^^) _ ^) e El, 

(4.2) - log^ V2nV6(A„(n-V3^2) - ^2^) G ^2 
= logP(A„(n-i/Vi) G A„(n-V3r2) G E2) 
= G„(n"^/^r; ai, . . . , a2r;bi, 62s) 

for the disjoint union of intervals 

r s 

(4.3) Ei:=[J[a2i-i,a2i] and E2 := [J[b2i-i,b2i], 

i=l i=l 
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with 



(4.4) ai = \/2^+ and bi = V2^+ 



The method here is to do asymptotics on (3.3) for n large. In the notation 
(1.14), define 

(4.5) L:=L,, + L^, E:=Eu + E^, 

with the understanding that t in £' now gets replaced by r. Setting k := n^^^ 
and changing variables, 

AiGn{n~^^^T; ai, . . . , a2r;bi,. . . , 62s) 

= AiHnir; -k{2e - V2ai), -k{2k^ - V2bi), . . . ) 

(4.6) 



fe\/2f ^ — + e"^/'' X! ^ • • ■:'^2r;Vi, . . .,V2r) 

All^^^Hn{T;Ui,. . . ,U2r;Vi,. . . ,V2r), 



b — >-v 



where the operators Ai and Bi are now expressed in u, v, r-coordinates, using 
the change of coordinates (4.4) and the chain rule. In these new coordinates, 
the Ai and Bi, Taylor expanded in 1/k, for large k, read as follows: 



)l, + o( ^ 



(4.7) 



^^ = ^'i''-{i-w^wy--yks 

A2 = 2k^ {l - ^L, + -L(^ _ 1 + ^r^L,) 
B2 = 2k* {l - + - 1 + 



12 
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Hence, from (4.7), 
^ B2A1 



L'-—{L + Lu)L 



1 



(4.8) 



+ ^{L{E - 2) + r2(4L„(L + L,) + LL,)) 
+ 



(^L{Eu - 2) + ^L^(^ + 2) + y (8LL„ + 18L„L^ + LL„)) 



1 



2A:2 



A:2 



fcn2 



11(1. 



l2 + L,L, + O 



1 



As will be shown in Proposition 7.1, using also (7.13), we have for u = 
(ui, . . . , U2r) and w = (ui, . . . , t;2s) (also for their derivatives with regard to 
the endpoints of the intervals) 



(4.9) 
with 



Hn{T; u, v)\n=ke = G{t; u, v) + o{l/k) for A; — > 00, 



G{t,u,v) :=logP(A(n) G^i,A(r2) ££2), 



which will be used below. Proposition 7.1 actually implies that the error in 
(4.9) has order A;~^logA;. Equation (3.3), with (4.7) and (4.8) substituted, is 
a series in k'^, for large k, but where the three leading coefficients, namely 
the ones of k'^, k'^ and k^, vanish: 



0: 



Ai/{V2k) 



(4.10) 



2r 

fc2 

2r 



/(V2k) 



{{Liu + L^'){Ly^Ey — L^Eu) +t'^{Lu — L^^Ly^L^^Hn 



1 



{{Lu + L^){LuEi, — LyEu) + t'^{Lu — Lt,)L„L^)G 
~ 2^{'^u ~ Lv)G, {Lu + Lv)^G} ^^^^^ 



using (4.9). In this calculation, we used the linearity of the Wronskian 
{X,Y}z in the three arguments and the following commutation relations: 



(4.11) [L„,KJ=L„, [Lu,E^,] = [Lu,L^] = [Lu,t]=0, [Eu,t] = t, 
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including their dual relations hy v; also {L^G, 1} = {L{Lu — Lv)G, 1} l 
It is also useful to note in (4.10), that the two Wronskians in the first ex- 
pression are dual to each other hy u<^ v. The point of the computation is to 
preserve the Wronskian structure up to the end. This proves Theorem 1.2. 
□ 

Proof of Corollary 1.3. Equation (1.15) for the probability 
G{t;u,v) :=logP(A(ri) < u,A{t2) < v), t = T2 -n, 
takes on the explicit form 

/ d'^ d^G r^d^G ^ d^G d'^G , ^ ^ 



dr V du'^ dv"^ J dv? dv V dv"^ ' du dv dv? ~'~ 

d^G f^d'^G d'^G d^G ^ 



^ ^ dv"^ du \ du^ ~^ du dv dv^ ~^ 

_^f^d__^d_\/d_ d_\ 

\ du^ dv dv^ du J \ du dv J 

This equation enjoys an obvious u ^ v duality. Finally the change of vari- 
ables in the statement of Corollary 1.3 leads to (1.16). □ 

5. The joint distribution for the Sine process. 

Proof of Theorem 1.4. Consider as in (1.13), the disjoint union of (in 
this case) compact intervals Ei := Ui=i[^2i-i5 ^2i] and E2 := [ji=i[v2i~ii 'V2i] ^ 
M. Then, again setting r = T2 — ti, 

Hn{T;Ui, . . .,U2r;Vi, . . . ,V2s) 

\ogpLl\ 2V^A(„,/2)+i (^^) G El, all 2V^A(„/2)+i G 



^^■^^ =logp(an A(„/2)+.(^) G^f,anA(„/2)+.(^) Gi?2^ 

= Gni —;ai, . . . , 02r; bi,. . . , b2s 

for the disjoint union of intervals 

r s 

(5.2) ^1 := |J[a2i-i,a.2i] and £'2 := |J [62^-1, &2i], 

1=1 i=l 

with 
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Note here G„ refers to the second formula in (3.4). Setting k :=n^/^, using 
the change of variables (5.3) and the chain rule, 

AiGniT/k'^;ai, . . .,a2r]bi, . . -Ms) 

= AiHn{T]2kai,...-2khi,...) 

= 2fc 1^ ^ — + e"^/'' ^ — j F„(t; Ui, . . . , U2r]Vi, V2r) 
= Al\^^^ HniT;Ul, . . . ,U2r;Vl, . . . ,V2r). 

b — >v 

In these new n, f , r-coordinates, the operators Ai and Bi, Taylor expanded 
in l//c for large k, read as follows: 

Ai = 2k(L- + ^ 

Bi=2k(L-^Lu + 0[ 

A2 = E-l-'^^{E,,-l) + 0[^^^^ 

B2 = i?-l-|(i?«-l)+o(l^^ 

Moreover, as will be shown in Proposition 7.3 (Section 7.3) we have the 
following asymptotic estimate (and similarly for its derivatives with regard 
to the endpoints of the intervals): 

(5.5) Hn{T\u,v)\n=k'i = G{t;u,v) + o{l) for/c^oo, 
with 

(5.6) G{t, u, v) := log P fall V27rSi f ^) G Ef, ah V2TrSi G E^) . 



(5.4) 




Using the expansions (5.4) of the Ai and Bi and later the commutation 
relations (4.11), yields the following for the Wronskian: 



1 



2, 



i2ky 



-{B2AiHn,BiAiHn + 2k'e--^' 

^2k((E-l-^iE^-l)](L- ^^L^ +o(^,]]H., 



{2k)H \\ k'^' vV k^y'w 



71) 



(2fc)' ((L-^L, UL-^L^ +Ol-,]]H^ 



2fc(L-(r/fc2)L„+0(l/fc4)) 
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-■^|(i?-l)LJf„,L2i/„ + i|^ 



Hence, subtracting the previous formula from its dual and using (5.5), 
1 



!^{E-l)LHn,L^Hn + ^^^ 



- |(2(^„ - 1)L + {E + l)L,)Hn, LHn + 
+ |(2(^, - 1)L + + l)Lu)Hn, Lif„ + + O (^^' 

{(£;-i)lg,l2g' + 1| 

- I (2(E„ - 1)L + + 1)L„)G, LG + i 
+ I (2(E, - 1)L + + 1)L„)G, LG + + o(l) 



Lu — Lv 
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Upon division by (L^G + ^)^, one finds 



'VL20+ 1/2/ 



+ {L'a + Ly] 



{2{E, - Eu)L + {E + l){Lu - L,))G 



L'^G + 1/2 

{{E - l)L + 2{E, - Eu)L + {E + 1)(L, - L,))G 



L'^G + 1/2 

-{E - l)L + 2{E, - Eu)L + iE + - L,))G 



2L, 



L2G + 1/2 
{2EyLu + {Ey — Eu — l)Ly)G 

iJgTT/2 

{2EuLy + {Eu — Ey — l)Lu)G 



2Ltj 

L-^G + 1/2 
In view of (5.6), the function 

G{t,u,v) = logP(all Siin) £ El, all ^^(ra) G E'^) 

= G(^^,V2TTu,V27rv 

satisfies (1.17) of Theorem 1.4. □ 

Proof of Corollary 1.5. Setting 

Ui=Xi + X2, U2 = Xi-X2, 

vi = yi+y2, V2 = yi-y2, 

the function 

(5.7) H{t;xi,x2;yi,y2) ■■=G{t;xi +X2,xi - X2;yi + y2,yi - y2) 
satisfies (1.18), ending the proof of Corollary 1.5. □ 

6. Large time asymptotics for the Airy process. This section aims at 
proving Theorem 1.6, for which we need the fohowing lemma: 



Lemma 6.1. The following ratio of probabilities admits the asymptotic 
expansion for large t > in terms of functions fi{u,v), symmetric in u and 
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v: 

P{A{^)<u,A{t)<v) _ ^ Mu,v) 
P{AiO)<u)P{A{t)<v) ' 

from which it follows that 

lim P(A(0) < u,A{t) <v) = P{A{0) < u)P{A{t) <v) = F2{u)F2{v)] 

this means the Airy process decouples at oo. 

Proof. This will be done in Section 7.5, using the extended Airy ker- 
nel. Note, since the probabilities in (6.1) are symmetric in u and v, the 
coefficients fi are symmetric as well. The last equality in the formula above 
follows from stationarity. □ 

Conjecture. The coefficients fi{u,v) have the property 

(6.2) lim fi(u, v) = for fixed veR 

u — >oo 

and 

(6.3) lim fi(—z, z + x) = for fixed x G M. 

z — ^oo 

The justification for this plausible conjecture will now follow: First, consid- 
ering the following conditional probability: 

P{A{0)<u,A{t)<v) 



P{A{t) < V I ^(0) < u) 



P{A{0)<u) 

fi{u,v) 



^2(^^) 1 + E t 

\ i>l 



and letting v ^ oo, we have automatically 

fi{u,v) 



1 = lim P{A{t) < V I ^(0) <u)= lim 



\ i>l ^ 



i>l 

which would imply, assuming the interchange of the limit and the summation 
is valid, 

(6.4) lim /,(n,?;)=0 

and, by symmetry 

lim fi{u,v) =0. 
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To deal with (6.3), we assume the following nonexplosion condition, whose 
plausibility is discussed in the Appendix: for any fixed t > 0, x G M, the 
conditional probability satisfies 

(6.5) lim P(A(t)>x + z\A(0)<-z) = 0. 

z — ^oo 

Hence, the conditional probability satisfies, upon setting 

V = Z + X, u = —z, 
and using lim^^^oo F2{z + x) = 1, the following: 

fi{-z,z + x) 



1 = lim PiA{t) <z + x\ A{0) <-z) = l+ lim ^ 

which, assuming the validity of the same interchange, implies that 

lim fi(—z, z + x) =0 for all z > 1. 

2^00 

Proof of Theorem 1.6. Putting the log of the expansion (6.1) 
G{t;u,v) = log P{A{0) < u,A{t)<v) 

(6.6) i>i ^ 
= logF2iu) + logF2{v) 

, fiiu,v) h{u,v)-ff{u,v)/2 

in (4.12), leads to: 

(i) a leading term of order t, given by 

(6.7) £^1 = 0, 
where 

\du dv J du dv 
The most general solution to (6.7) is given by 

hi{u,v) =ri{u) +r3{v) +r2{u + v), 
with arbitrary functions ri,r2,r3. Hence, 

P(^(0) < n, A{t) <v) = F2{u)F2{v) {l + + . . >j 

with hi{u,v) = fi{u,v) as in (6.1). Applying (6.2), 

ri(u) + r3(oo) + r2(oo) = for all ti G M, 
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implying 

ri{u) = constant = ri(oo), 

and similarly 

r^(u) = constant = r^(oo). 

Therefore, without loss of generality, we may absorb the constants ri (oo) and r3(cxD) 
in the definition of r2{u + v). Hence, from (6.6), 

fi{u,v) = hi{u,v) = r-2{u + v) 

using (6.3), 

= lim fi(—z,z + x) = r2 (x) 

z — >oo 

implying that the /ii(u,f)-term in the series (6.6) vanishes. 

(ii) One computes that the term /i2(ti, v) in the expansion (6.6) of G{t; u, v) 
satisfies 

(6-9) Ch, = —^—^ - —^—^ with g{u) ■.= \ogF,{u). 

This is the term of order , by putting the series (6.6) in (4.12). The most 
general solution to (6.9) is 

h2{u,v) = g'{u)g'{v) +ri{u) +r^{v) +r2{u + v). 

Then 

P(^(0)<n,A(t)<t;)=e^(*'"''') 

hi{u,v) 



F2{u)F2{v) expy^ 



^ , ^ , N / h2{u, v) \ 
F2{u)F2{v)^l + ^^ + --j. 



In view of the explicit formula for the distribution F2 and the behavior (1.6) 
of q{a) for a /" 00, we have that 



lim g'(u) = lim (\ogF2(u))' 

= lim / q^(a)da = 0. 



Hence 

0= lim f2{u,v)= lim h2{u,v) = ri{oo) + rs{v) + r2{oo), 
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showing r3 and similarly ri are constants. Therefore, by absorbing ri{oo) 
and r3(oo) into r2{u + v), we have 

/2(n, v) = /i2(n, v) = g'{u)g'{v) +r2{u + v). 

Again, by the behavior of q{x) at +00 and —00, we have for large z > 0, 



00 roo 



g'{-z)g'{z + x)= q\a)da q\a) da < cz^/^e-^'/^ 

J—z J z-\-x 

Hence 

0= lim /2(— z, z + = r2(x) 

z — ^00 

and so 

f2{u,v) =h2{u,v) =g'{u)g'{v), 

yielding the term in the series (6.6). 

(iii) Next, setting 

G{t;u,v) = log P{A{0) < u,A{t)<v) 
(6-10) , , , , ^ g'iu)g'{v) h,iu,v) 

= 9{u) + 9{v) + ^2 + — p + ■ ■ ■ 

in (4.12), we find for the term: 

Chs = 0. 

As in (6.7), its most general solution is given by 

h3{u,v) =ri(n) +r3{v) +r2{u + v). 
By exponentiation of (6.6), we find 
P{A{0) <u,A{t) < v) 

= F2{u)F2{v) (1 + + p + • 

The precise same arguments lead to h-3{u,v) = 0. 

(iv) So, at the next stage, we have, remembering g{u) = logi*2(ii) 

(6.11) G{t] u, v) = g{u) + g{v) + ^ + — — + • • • 

with 

(6.12) fi{u, v) = hi{u, v) + \hl{u, v) = h^iu, v) + ^g' (u)"^ g' {vf . 
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Setting the series (6.11) in (4.12), we find for the term: 

Ch -2(—(—Y d^d^fdi dg\ 

^ \ du^ \ dv^ J dv^ \ du^ ) ) du^ dv^ \du dv ) 

2 V du"^ dv\dv J dv'^ du\du) J 

(6-13) \du^ dv"^ dv^ du^ ) \du^ dv dv^ du ) 

= 2{2q{u)q'{u){q{v)q'{v) + l) 

roc 

-q{u)q"{u)q\v)-{q'{u)fq\v)) / q^ 



+ 2q{u)iqiu)q'{v)q"iv) + q iu)q{v)q" (v) - 2q{u)q\v)q' {v)) 
— same with v. 

The latter is an expression in q{u), q{v) and its derivatives and in Q^{(y) da 
and /^°° (f'ioi) da. It is obtained by substituting in the previous expression 



g{u) = I {u — a)q {a) da 

J u 

and the Painleve II differential equation for q{u), 

uq{u) =q{u) -2q{uf, 
in order to eliminate the explicit appearance of u and v. Now introduce 



giu) = 


poo 

1 {u — a)q'^{a) da, 

J u 


9i{u) = 


roo 

/ {u — a)q'^{a) da, 

J u 


92{u) = 


roo 

/ {u — a)q'^{a) da. 

J u 



Note g'{u)=j:^q\a)da,g"{u) = -q\u) and g[{u) = q'\a) da,g'^{u) = 
q'^{a) da. The most general solution to (6.13) is given, modulo the null- 
space of £, by 

hAiu,v) = ^{g" {u)g' {vf + g" {v)g' {uf + g" {u)g" (v)) 

+ g'(u){2giv)+g[iv)-g'^iv)) 
+ g'{v){2g{u)+g[{u)-g'^{u)) 

+ / q'^{a)da / {2{v — a)q'^{a) + q'^{a) — q'^{a)) da 



+ same with n <-> f . 
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This form, together with (6.12), imphes for the function fi{u,v): 

fi{u,v) = hi{u,v) + \g' {uf g {vf + ri{u) + r:i{v) + r2{u + v) 
= ^ai{u)bi{v) +ri{u) + r^iv) +r2{u + v). 

i 

Using the asymptotics of q{u), one finds 

I CG ^ U — ^ OC 

ai{u),bi{u) < <^ I ,3' ' 
[ c\u\ , u —00, 

and so, by the same argument, 

ri(u) = r2{u) = r3(u) = 0. 

Therefore, we have 

U{u, v) = hi{u, v) + \g{ufg{vf, 

with h/^{u,v) as in (6.14), thus yielding (1.19). 

Finally, to prove (1.20), we compute from (1.19), after integration by 
parts, taking into account the boundary terms, using (1.6), 

E{A{{))A{t)) = f f uv-^P(Am < u, A(t) < v) du dv 
J Jm? ou ov 



uF^{u)du / vF^{v)dv 

D J — 00 

roc poo 

+ -2/ F!,{u)du F^{v)dv 

r J~oo J— CO 

+ -I [ [ {^{u,v) + ^{v,u))dudv 

= (£^(^(o)))^ + ^ + ^ + o(^), 

where 

c:= / / {^{u,v) +^{v,u))dudv = 2 / / ^{u,v)dudv, 

J JR2 J J^2 

thus ending the proof of Theorem 1.6. □ 



7. The extended kernels. The joint probabilities for the Dyson, Airy and 
Sine processes can also be expressed in terms of the Fredholm determinant 
of matrix kernels, the so-called extended Hermite, Airy and Sine kernels 
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(considered in [6], [11] and especially in [12] and [8]), defined for subsets 

Ei cR, 

with Kt^tj being one of the following kernels: 

oo 

E ^~''^*'~*'Vn-fc(2;)l/3n-fc(?/), if ti>tj, 

L 



e''^^^-^''^(Pn-k{x)(Pn-k{y), iiti<tj, 



fc=i 



A:=— oo 

oo 



-2{ti-ii)Ai(x + z)A\{y + z) dz, 



if ti > tj, 



e^fe *i^Ai{x + z)Ai{y + z)dz, if ti<tj. 



1 r 



^^y^ cos z{x - y)dz. 



if ti > tj. 



^-z^itj~U)/2 



where 
(7.2) ^k{x) 



e ^^^'^Pk{x), for A; > 0, with pk{x) 



cos z{x — y) dz, if tj < tj , 
Hk{x) 



0, 



for A; < 0; 



2^-/2VM7ri/4' 



Pk{x) are the normalized Hermite polynomials, and Ai(x) is the Airy func- 
tion. Now we make a few comments about these kernels. 

7.1. The Fredholm determinant of extended kernels. Letting x{t) denote 
either the largest eigenvalue A„(t) in the Dyson process, or the Airy process 
A(t), or the collection of S'j(t)'s in the Sine process, the probability is now 
defined by (drop the superscripts in K^'"',K'^,K^) 

P{x{ti) e Ei,l<i<m) 

= det(/ - z{Kt^tj)i<i,j<m)\;,=i 

0<ri<N ^ I 1 
1 r, = iV 



1+E 



da 



(m) 



N=l 



X det ( {Kt^ti (af ^ , af)) i<,<, „ , ^ ^ , 



z=l 
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where the A^-fold integral above is taken over the range 



-OO < a^-^^ < • • • < Ori^ < oo 



-oo < < • • • < drj < oo _ 



with integrand equal to the determinant of an x matrix, with blocks 

y, ■ 



given by the x matrices {Kt^ti{o(f'\ af^)) i<i<rf, . In particular, for m 



2, we have 

Pix{ti)eEi,x{t2)£E2) 



N=l 



X det 



0<r, a<N 
r+s = N 



{ 



/ 1 



Y[ dai J]^ dj3i 



l<i,j<r 



{Kt^t2{oii,l3j))i<.< 



l<j<a 



l<j<r 



Z = l 



7.2. The extended Hermite kernel tends to the extended Airy kernel. Given 
the substitution 



u 



(7.3) Si 



we have: 



t 



X I— > \/2n + 1 H ^ — — 



n 



1/3'' 



n 



1/3' 



V2n + 1 + 



^/2nV6 



Proposition 7.1. T/ie extended Hermite kernel tends to the extended 
Airy kernel, when n oo, uniformly for u,v (z compact subsets C M: 

\Ktf{x,y)dy\g^ - K^^,{u,v)dv\=o(^^^^. 

Before proving this proposition, we need the following estimate: 

Lemma 7.2. For large n > and — Mon^/^logn < k < Mgn^/'^ logn, with 
fixed Mq > 0, we have 



(7.4) <^„_J V2^^ + 



2i/4 



n 



k 



73 



V2ni/6^ 

wii/i i/ie following uniform bound in compact subsets M.: 

' log n ' 



\Eu{k,n)\<0 



n 



2/3 
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1-6 1+5 
Fig. 1. 

Proof. Here one needs the asymptotics for the Hermite polynomials 
when z G 6*2,5, as in Figure 1; it is given by (see [3]): 



(pn{zV2n) =pnizV2n)e ""^ 
(7.5) = (2n)V4 Mfniz)) 



z-l 



fniz)) Ai'(/„(Z)) 



the error term being uniform in C2,5q for some 6q > 0. This captures the case 
> in the statement of Lemma 7.2. The case u <0 would be captured by 
a similar estimate valid in the region Ci^s- To explain (7.5), the equilibrium 
measure for the Gaussian distribution is given by the well-known Wigner 
semicircle 

TT 

and 

2/3 



/„(z) = -(-n)2/3^ 

{zv 1 — z" — arccos z 
Setting z = 1 + X for small x>0, one computes 

/„(l + x) = 2n2/3x(^l + ^ + ---) 



and 



fn[z) 



z=l+x 



V2ni/6fl + Ax + -. 



.z-l 
Defining x such that 

V2^+ = ^2{n - k){l + x). 
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one computes, for k = Mv}^^ logn and \M\ < Mq 
(7.6) 



1? 



u 2k + 1 ku 



loff n + u 1 Mu log n 



2n2/3 4n 4^4/3 

Thus for X behaving as (7.6) and for k = Mn^^^Xogn, we deduce from the 
formulae above 



fn^k{z)l=i+. = 2(n - kf^x ( 1 + + • ■ 



2k 



J—^fn-k{z) 



z=l+x 



V 3n A 10 
V2{n-k)^/'^(^l + ^x + -- 

y2„■/6('l+o('^^i5|!! 



and 



(n - A;)V4 ni/4 



-, ^ 
1 + — + 

4n 



n 



1/4 



l + O 



Mlogn 

2/3 



n 



Using the asymptotics (7.5), one computes for k = Mn^/^ log n and |M| < 
Mo, 



fe V2n + 1 + 



V?„_fc(y^2(n-/c)(l + x)) 

l + Ojl/jn-k)) f fz + l 
(2(n-A;))V4 \U-1 

z + 1 



l + 0(l/n) 



1 + 



(2n)V4 



z- 1 
Mlogn 

^2/3 

Af log n 

^2/3 



fn~k{z) 



1/4 



Ai(/n-fc(;j)) 
-1/4 

Ai'(/„-fc(2)) 



2=1+0; 
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X Ai ( M log n + u + 
l + O 



1 



1 



\/2nV6 



Mlogn 



+ 



M log n 

^2/3 



2/3 



X Ai'(Mlogn + M + 



1 



2nV3 



+ 



logn 
7^ 



21/^(1 + 0(Mlogn/n2/3)) 



nVi2 
M log n 

2/3 



A-/,x, N Ai'(Mlogn + u, ^ 
Ai(Mlogn + u) H ^ „ f,, ^ + 



Mlogn 

2/3 



n 



2nV3 
Ai'(Mlogn + n) 



2i/4 

2i/4 
n 



1/12' 



Ai{u + Mlogn 
Ai u + 



1/3 



n 



l + O 



l + O 



2nV3 
M log n 

^2/3 

M log n 



+ 



Mlogn 

j^2/3 

M log n 

^2/3 



n 



2/3 



ending the proof of Lemma 7.2. □ 



Proof of Proposition 7.1. As a first step, in a recent paper Krasikov 
[10] shows the following inequality for A; > 6 and for a universal constant c: 

1/2 

(7.7) max\Hk{x)\e 



/ 2W4k-2 



\k/2)\ VfcV6V8A;2-8A; + 3. 
From Stirling's formula n! = ^/2^m n"e~"(l + 0(l/n)), it follows that 



2'=/2(A;/2)! 



2 \V4 
irk 



1 + 



This estimate, estimate (7.7) and formula (7.2) show that, for k > some 
fixed kf) and some constant c'. 



(7.8) 



1 



max |vp,(x)| = ^— ^=-^ max |/ffc(x)|e 



-x2/2 < 



A;i/i2' 



Using both estimates (7.4) and (7.8) in (7.1), one computes for t> s, taking 
into account substitution Si, as in (7.3), 



H,n, 
t,s \ 
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[Mn^/S log n 



k \/ 21/4x2 dv 



k=0 



xAi[v + 



< 



ni/3;VnVi2; ^/2ni/6 



k=0 



X 9J„-A: V2n + 1 + 



[Mn'^/^ log n 



(7.9) 



2i/4 



fe=0 



[Mn^/^ logn] 
+ ^ g-(fc/nV3)(j_,) 

k=0 



nVi2 

21/4 



Ai(t; + -^)(l + ^,(fe,n)) 



n 



1/6 



Ai I ti + ) Ai ( V + 



n 



1/3 



nV3 



dv 



X |l-(l + K(A:,n))(l + ^^(A;,n))| 



< 



E 



-(A;/ni/3)(t-s) 



fc=[Mnl/3log n]+l 



X (/?„_fc V2n + 1 + 



U 



X V5„_fc V2n + 1 + 



dv 



< 



^/2ni/6 



n-l-[Mn^^^\ogn] 
J2 g-(n-m)(t-s)/ni/3 _l_ 

mV6 



E 

m=0 



+ c"Mfi^(l+o(l); 



where c" is determined by the maximum of the Airy function Ai(2;) on the 
semi-infinite interval (0, cxo). 
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Setting n' = n — 1 — [Mn^^^ log n] , the sum in the last expression is esti- 
mated as follows: 



1 1 



76 



-{n-m){t-s)/n^/^ 



1 1 " ' 

< ^ ^ -{n-m){t~s)/n^ri i ^ {n-m){t-s)/n^/'i 



1 1 



-{n-£+l)(t-s)/ni/3 _ ^_(„+i)(t_s)/„l/3 



_^£-l/6(^g-(n-n')(t-s)/nl/3 _ ^_(„_£+i)(t_s)/ni/3 ^ 

1 /3 

_ ^ r„^(n-^+l)(f~.)/nV3 _ („+i)(t_^)/„i/3 

Picking £ = 0(n/2), all terms above tend to exponentially fast, except the 
term (*), which requires some attention. Choosing n' = n— [Mn^/^ logn] — 1, 
so that n — n' = 0(n-^/'^ log n), the coefficient of that term is 0(n~^^(*~*)). 
Therefore that term gets small, when n — > oo and £ = 0(n/2), 

1 ni/3 

.O(n-*^(*-))^0. 



^/2ni/6 ^1/6 



The proof is ended by observing that the second term in the first difference 
of (7.9) is a Riemann sum converging to the extended Airy kernel, that is. 



[Mtt}/^ log n] 

E =-"='"""'"-"Ai(»+;;?s)Ai(. + ;^ 

fe=0 



i/3vf-.^../ K \ ,./ k \/2^/'^\^ dv 



y nyy\n^/i2j ^^1/6 

= dv e""(*-^) Ai(n + z)Ai{v + z)dz + ( ^^^"f' ^ ■ 

This establishes the convergence for t> s. 

For i < s, one computes, again using in (7.1) the estimates (7.4) of Lemma 
7.2 and (7.8), 

K^s''{x,y)dy^^^ 

[Monl/^logn] , r \ / r \ / ol/4 \ 2 j 

V e-(^/"''')(^-*)Ai(^n--^^Ai(^;;--^V^^ 

k=[j * 
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< 



OO / \ 



X ipn+k i V2n + 1 + 



log n 

E e 

k=0 



-{k/n^/'^)(s-t) 



2l/4 
77712 

2I/4 



Ai u 



A; 



fc=0 



nVi2 
2I/4 



Ai u 



ni/3 
A; 



n 



1/3 



(l + £;(A:,n)) 



{l + E{k,n)) 



dv 



n 



1/6 



(7.10) 



X \l-{l + Eu{k,n)){l + E^{k,n))\ 



ni/3 



^/2ni/6 
dv 



^/2ni/6 



< 



+ E 



fc=[Mni/3log n]+l 



-{fc/ni''3)(s-t) 



iPn+k [ V2n+ 1 + 



n 



< 



X V3„,4.fc ( V2n+ 1 + 



E 



V2ni/6 



V2ni/6 



V2ni/6 



^ \ fc=[AfnV3logn]+l ^ ^ ^ 



dv 



< I g/2^-M(s-i)-l/6 ^-(m/nl/3)(,_i) 



E^ 

m=l 



The rest of the proof goes the same way as before, endmg the proof of 
Proposition 7.1. □ 
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7.3. The extended Hermite kernel tends to the extended Sine kernel. Given 
the substitution 



(7.11) 
we have: 



1 1 



, S I— > ,X I 

2n 2n 



iru TTV 



2n 



Proposition 7.3. The extended Hermite kernel tends to the extended 
Sine kernel, when n — > oo, uniformly for u,v G compact subsets C M: 



S2 



Sketch of proof. From [13], page 198, it follows that for |x| < M, 



1 



2'=/2(A;/2)!7ri/4 



(-^] ^ (cos( V2k + Ix - ^] + O (-^ 



cosf x\j2k + 1 



k-K 



+ 



6V2A; + 1 



sinl x\j2k + 1 



/cvr 



\k-K'^ ) 



\Vk 



Using the substitution ^2 as in (7.11), one computes in (7.1) for t>s, 



K^f{x,y)dy\, 



-7r2(t-s)/2_ 



n-1 

7r2(t-s)fe/2n 



(7.12) 



-7r2(t-s)/2^ 



7r« 



E 

fc=0 



2n 
vrn 



■KV 



'2n 



2n 

■KV 

y/2n 



dv 



dv 



TT'^{t-s)k/2n 



X e 
(I) + (11) 



y^k 



E+ E 

fc=0 k=a+l 
TTU 



/2n 



/2n 



(if 



where a > is the minimal integer above which Krasikov's estimate (7.8) 
holds. But then expression (I) tends to exponentiahy, when n tends to oo, 
and (II) is estimated as follows: 



(11) < 



/2n , A;V6 



2n 
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which tends to for £,n ^ oo such that tl^jn^l'^ 0. Then the sum ap- 
pearing at the first hne of (7.12) can be estimated as fohows: 



'2n 



fc=£+l 
n-1 



vrn \ / irv 

ipk I dv 



-■K^(t-s)/2 



E 

k=e+i 



^Tr'^k{t-s)k/{2n) 



n 



n 



1/2 



/ k 1 nvr k 

cos TTU\ — i 

I y n 2n 2 n 



/ k 1 mTk\ ^/l 

X cos -KV\ 1 +0 \ —pz 

\ \ n 2n 2 n 



n—l / 
-7r2(t™s)/2 J_ ^TT^ k{t- s) / {2n) 



k=i+l 



2n 



\k 



1/2 



COS iriu — v)\ — I 

\ \ n 2n 



. , k I k' 

+ COS 7r(M + f)W — I nvr— 

V \ n 2n n 



dv 



^ _^g7r2(x-l)(t-s)/2 
'i/n 2yfx 

cos(7r(ii — f ) 

+ cos(7r(n + v)\fx — nirx) + O + ^ ( 



dv 



dv 



1 r 



cos(z(n — w)) + cos(z(u + t>) — nz^ /tt) 



1 
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TT Jo 



when I and n tend to oo. The integral involving the second cosine in the ex- 
pression (*) above is an oscillatory integral and thus tends to zero faster than 
any power of n by the Riemann-Lebesgue lemma. The case t<s proceeds 
along similar lines, establishing Proposition 7.3. □ 

7.4. Convergence of Fredholm determinants and their derivatives. Here 
we give a schematic argument, based on the customary formula logdet(I — 
K) = Trlog(/ — K), where X is a kernel restricted to a disjoint union of 
intervals E and where K + 6K tends to K. Then, from 

, . det(/ -K-5K) = det(/ - K) det(I - (I - KY^5K) 

^ ' = det(/ -K){1- TV(/ - K)~^6K + o{5K)), 

one sees that det(I — K — SK) tends to det(/ — K), when 6K tends to 0. 
Also, given pi, ■ ■ ■ ,P2r the endpoints of the set E (see [7, 15]), 

-^logdet(/ -K) = {-lf-^K{I - Kr\pk,Pk), 

OPk 

where here "=" means "kernel of," evaluated at {pk,Pk) and so 
■^logdet{I-{K + 6K)) 

OPk 

= {-1)''-\K + 5K){I -K- 5K)-\pk,Pk) 
= {-lf~\K{I-Kr\pk.Pk) 
+ {5K{I-K)-^ 

+ K{I - K)-^5K{I - K)-^){pk,Pk) + o{SK)] 

d 

= —logdeiil - K) + 0{6K). 

OPk 

Since by Propositions 7.1 and 7.3 the extended Hermite kernel converges to 
the extended Airy and Sine kernels, this argument shows the convergence 
of the corresponding Fredholm determinants and their first derivatives with 
respect to the end points of E. In a similar fashion one proves the result for 
higher derivatives. 



7.5. An a priori asymptotic expansion for the joint Airy probability. The 
proof of Theorem 1.6 in Section 6 was based on an a priori asymptotic 
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expansion for the ratio below in 1/t for large t = t2 — ti. This can be found 
in [18] and proceeds as follows: 

P{A{ti) < u,A{t2) <v) _ det(/ - {Kf^t,)i<i,3<'^) 
P{A{ti) < u)P{A{t2) <v)~ det(/ - i-^^tj det(/ - K^) 

(T.14) M'-iZ "o")) 

= det(/ - /C12/C21) 

_-, I fi{u,v) 
fi ' 
i>l ^ 

where 

/C12 := XKoo)(2;)i<'TO(2;,y)X[«,oo)(?/))~"^X[M,oo)(a;)-f^^t(x,y)x[,;,oo)(y), 

/C21 := {I - X[v,oo)(.x)KQQ{x,y)x[v,oo){y)y^X[v,oo){x)K^o{x,y)x[u,oo){y), 
with 

K^o{x,y)= / Ai{x + z)Ai{y + z)dz, 
Jo 

K^q{x, y) = / e-"*Ai(x + z)k\{y + z)dz = 0{l/t), 
Jo 

rO 

K^^{x,y) = - / e'^ Ai{x + z) Ai{y + z)dz = {I /t). 
J— 00 

APPENDIX: REMARK ABOUT THE "NONEXPLOSION" 
CONJECTURE 

To discuss the conjecture (1.21), consider the Dyson Brownian motion 
(Ai(t), . . . , Xn{t)) and the corresponding Ornstein-Uhlenbeck process on the 
matrix B. Then, using the change of variables 

M,; : 



V(l-c2)/2' 

and further M2 1— > M : = M2 — cMi in the M2-integrals below and noting that 
max(specMi) < —z and max(specM2) > a imply max(spec(M2 — cMi)) > 
a + cz, we have for the conditional probability, the following inequality: 

P(A„(t)>a| A„(0)<-z) 

r c-(l-c2)TrM2/2 f j/i^ ^- Tr(Af2-cMi)2/2 

_ Jmax(specMi)<-^ Q^lg ^ ^ Jmax(spec M2)>a "'^S e ^ " 

Jmax{spccA/i)<-2 "^'^^1 ^ ^ JM2e'W„ ^^'^S e ^ ^/ 

f ^A/f p-(l-c2)TrM2/2 f j)i^„-TrM2/2 
Jmax(specMi)<-2 "^'^-'1 ^ ^ Jmax(spec M)>a+cz "^''^ ^ 

Jmax(spccA'/i)<-2 "^'^1 ^ JM&Hn,"'^^^^ 

= P{\n{t)>a + cz), 



< 
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implying 

lim P(Xn{t) > a I An(0) < -z) = 0, 

z — ^oo 

and a fortiori, 

lim P{Xn{t) >x + z\ A„(0) < -z) = 0. 

z — >oo 

It is unclear why the limit (6.5) remains valid when n ^ oo, using the Airy 
scaling (1.4). 
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